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Abstract 

We describe a light-cone QCD sum rule calculation of the B s — > /o(980) transition form factors useful to predict the 
branching ratios of the rare decays B s — > f)( + £~, B s — > f)vv and of B s — > J/if/fo decay assuming factorization. We 
compare this channel to B s — > J/if/</> as far as the possibility to determine the B s mixing phase is concerned. 

Keywords: B s decays, QCD sum rules, CP violation 



1. Introduction 

Rare B s decays induced at loop level in the Stan- 
dard Model (SM) are sensitive to new physics (NP) 
effects that may enhance their small branching ratios 
(]]]. Besides, the analysis of the B s unitarity trian- 
gle of Cabibbo-Kobayashi-Maskawa (CKM) elements: 
VusV* b +V cs V* b +Vt s V* b = provides an important test of 
the SM description of CP violation. One of its angles, 
Vi v* 1 

B s = Arg - y yt , is expected to be tiny in the SM: 

B s ^0.019 rad. Recently CDF H and DO H Collabo- 
rations have indicated larger values with sizable uncer- 
tainties, although the latest CDF analysis |4J] seems to 
reconcile the SM with data. Hence the precise measure- 
ment of B s is a priority for forthcoming experiments. 
In this paper we describe the light-cone QCD sum rule 
(LCSR) calculation of the B s / (980)Q form factors 
||5[1, using the results to predict the branching ratios of 
the decays B s -> fa£ + {~, B s -> f vv in the SM. We also 
study the mode B s — > J /if/ fa that allows to access B s @]. 



2. B< 



/o form factors in Light-Cone Sum Rules 



The matrix elements involved in B s — > fo transitions 
can be parameterized in terms of form factors as 
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(MPf )\so- MV y5q v b\B s (p Bs )) = 
FAq 2 ) 

TUB, + »l/ c 

where P = pg s + p/ a and q = pg s -p/ - To compute such 
form factors using light-cone QCD sum rules (LCSR) 
10] we consider the correlation function: 



i J" d 4 x ( 



"(foMT^W, ^(0)110) (3) 



where jV, is one of the currents in the definitions (Q])- 
d2J of the form factors: jr, = Jf, - sy^y^b for F\ 
and F , and j F] = = scr^yrfb for F T . j Tl = 
biyss interpolates the B s meson; its matrix element be- 
tween the vacuum and B s defines the fs t decay constant: 

(B s (p Bs )\biy 5 s\0) = ^Jj/a,. The LCSR method con- 
sists in evaluating the correlator Q both at the hadronic 
level and in QCD. Equating the two representations 
gives a sum rule suitable to derive the form factors. 

The hadronic representation of the correlator in (01 
can be written as the contribution of the B s plus that of 
the higher resonances and the continuum of states h: 



n H (p /o ,<7) = 



(MPfJj^BApfy + q)){B s ( Pk + q)\M\0) 



m l s ~ (Pfo + <?) 2 
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•J Sn 



ds- 



p h (s,q 2 ) 
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'so ' s-ip/o + q) 2 ' 

where higher resonances and the continuum of states 
are described in terms of the spectral function p h (s, q 2 ), 
which contributes starting from a threshold so- 
To evaluate the correlator in QCD we write it as 

Imn QCD (>, g 2 ) 

* - (pa + q) 2 ' 

August 23, 2010 



n QCD (p/o,4) = - ds 

H J (m/j+m,) 2 



Expanding the T-product in (0 on the light-cone, we 
obtain a series of operators, ordered by increasing twist, 
the matrix elements of which between the vacuum and 
the /o are written in terms of fo light-cone distribution 
amplitues (LCDA). Since the function p h in © is un- 
known, we use global quark-hadron duality to identify 
p h with p QCD = ±ImIi QCD when integrated above s 10: 



f 
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i r° 
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Using duality, together with the equality FL n (pf ,q) 
n QCD (p f ,q), we obtain from Eqs. © and ©: 

(MPfJjrMiPfo + q))(Bs{p fo + q)\jrM 



If • 
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perform 
two sides 
1 



+ Q 2 )" 



a Borel transformation of 

in ©, exploiting the result 

exp(-s/M 2 ) 2 2 

where Q = -q 



(M 2 )" (n - 1)! 
and M 2 is the Borel parameter. This operation improves 
the convergence of the series in Tfi CD and for suitable 
values of M 2 enhances the contribution of the low lying 
states to I1 H . Applying it to n H and rfi CD we get 



(MPfMr, \Bs{pBMBs{p Bs )\ir 2 |0)exp 
1 f s ° 

ds exp[-s/M 2 ] Imn QCD (s,^ 2 ) 

f J(m b +m s ) 2 



M 2 
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(Pb s - Pf + q)- From O we derive the sum rules for 
Fi, Fq and Ft, choosing either yV, = Jf, or jr, = J^ T - 

In the calculation of n^ CD we consider fo as a ss 
state modified by some hadronic dressing [9]. Possi- 
ble fo — a mixing may only affect the overall nor- 
malization of the form factors at zero recoil, a system- 
atic uncertainty in our numerical results. We refer to 
J5I] for the definitions of the fo LCDA, for numerical 
values of the input parameters as well as for the final 
expressions of the form factors obtained from O. We 
fix so = (34 + 2) GeV 2 , which should correspond to the 
mass squared of the first radial excitation of B s . As for 
the Borel parameter, the form factors for each value of 
q 2 depend on it. The result is obtained requiring sta- 
bility against variations of M 2 . In Fig.Q~]we show the 
dependence of F\{q 2 - 0) on M 2 . We observe stability 
when M 2 > 6 GeV 2 , and we fix M 2 = (8 ± 2) GeV 2 . 

To describe the form factors in the whole kinemati- 
cally accessible q 2 region, we use the parameterization 



Table 1: Parameters of the B s — > fo form factors by LCSR. 
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' (q ^ 1 - ajq 2 /m 2 B ^ + bj(q 2 /m 2 B ) 2 ' 
collect in Table Q] the parameters Fi(Q), a, and b\ ob- 
tained fitting the form factors computed numerically. 
The q 2 dependence is shown in Fig. [2] The uncertainties 
in the results are due to the input parameters, so and M 2 . 
The parameters a, and bi are close for F\ and F j. The 
reason is the following. In the heavy-quark limit and in 
the large energy (LE) limit of the recoiled meson, the 
Bs — > fo form factors can be related [11] as follows: 



— F T (q 2 ) = Fdq 2 ) = ^F Q (q 2 ), (8) 
mg s + rrif a lb. 

withg 2 = m 2 B -ImsE (neglecting m 2 o )- The first equal- 
ity in (O predicts the same q 2 dependence for F\ and 
Ft in the LE limit. For the parameters of Fq, the second 
equality gives: ao = -1 +a\, bo — 1 —a\ + b\ which, us- 
ing the results for a\ and b\, gives a Q LE ^ 0.44 + 0.1 and 
b Q ^ 0.15 + 0.12. Hence the first relation is respected 
in our calculation, while not much can be said about the 
second one due to the uncertainty affecting bo- 



3. Semileptonic B s -* fot + { and B s -* / vv decays 

B s decays induced by the b — > s transition can con- 
strain new Physics scenarios. For example, they are 
sensitive to the compactification radius of universal ex- 
tra dimensions ifllh . Among such modes we consider 
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Figure 1: Dependence of F. ' °(0) on the Borel parameter M 2 . 
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Figure 2: dependence of the B s — > Jo form factors. 



— > /o^ + ^ and B s — > /ow, using the B s — > /o form 
factors to compute their branching ratios. 
The SM effective Hamiltonian describing the transition 
b -> str is: 



tial decay width of B s — > fa( + ( reads: 



dT(B s f ( + n G\ a 2 em \ V th \ 2 \V; s \ 2 q 2 - Am) VI 




5\2m\ 7T 5 



3q 2 



6m 2 (ml - m 2 f /F 2 (q 2 ) + (q 2 - 4m 2 )AF 2 (q 2 ) 



+(q 2 + 2m 2 e )A 



C 9 F x (q) + 



2Cf(m b -m s )F T (q 2 ) 



m B s + m fo 



with A = A(m 2 B ,m 2 fo ,q 2 ) = (m^ - q 2 - m 2 Q ) 2 - 4m 2 fo q 2 , 
a em the fine structure constant and mt the lepton mass. 
Analogously, the effective Hamiltonian for b — » svv is 



G F ae m VtbV? s 

V2 2n sm 2 (6 w ) 



rj x X{x t )0 L = C L L 



(10) 



where Ol = (sy^il - ys)b) (vj^(l - ys)vj and W is the 
Weinberg angle; the function X(x,) (x, - mf/m 2 , with 
m t the top mass and mw the W mass) has been computed 
in |Q3] and 00, while ^ 1 El 53, El- From 
Ht^svv the differential decay width is obtained: 



dF(B s -> / vv) 
dq 2 



\C L \ 2 A 3 ' 2 (m 2 B ,m 2 q 2 ) 
3 , ' . /0 IF!^ 2 )! 2 



96m B 7r 3 



Referring to [5] for the values of the parameters, we get: 

5^-1) x i(T 8 

I +0.4x ^ , n-8 



fofr) = (9.5!^) x 10- 



&R{B 

mBs -» for + T~) = (i. i :;;;!) x io 



(11) 



(8.7!^) x 10~ 7 



with I = e, \i. Hence these decays are accessible at the 
LHCb experiment at the CERN Large Hadron Collider 
and at a Super B factory operating at the T(55 ) peak. 
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= - 4 -^VtbK J] CAn)0&i) , (9) 



Gf being the Fermi constant and Vh the elements of 
the CKM mixing matrix (we neglect terms proportional 
to V u bV* s ). The expression of the operators (9, can be 
found e.g in 01311 . The Wilson coefficients in (0 are 
known at NNLO in the SM [14]. C 3 - C 6 are small, 
hence the contribution of only Oj, Og and 0\q can be 
kept for the description of the b — > s( + €~ transition. We 
use a modified Ct , which is a renormalization scheme 
independent combination of C7, C% and C2, given by a 
formula that can be found, e.g., in lfl5ll . 

The matrix elements of the operators in Hb^ s e+e- can 
be written in terms of form factors, so that the differen- 



4. Nonleptonic B s -» J/ifffo transition 

In the B s sector, B s — > J/4"f> is the golden mode to in- 
vestigate CP violation. Analysing it, CDF [2] and DO 
^ Collaborations have obtained values of the B s mix- 
ing phase <p s = -2B S much larger than expected in the 
SM, modulo a large experimental uncertainty. Hence, 
it is of prime importance to consider other processes to 
measure B s , as B s — > J/if/T]^ and 7/i/r/o(980) in which 
the final state is a CP eigenstate and no angular analysis 
is required to disentangle the various CP components, 
as needed for B s — > J/i(/(p. However, the reconstruction 
of B s modes into 77 and 77' is experimentally challenging, 
since the subsequent 77 or 77' decays involve photons in 
the final state. The case of fit seems feasible, since /o 
essentially decays to 7T + 7r~ and to 2tt° lfl9ll . 



3 



From the theory viewpoint, the quantitative descrip- 
tion of nonleptonic decays is challenging. Using the 
operator product expansion and renormalization group 
methods one can write an effective hamiltonian as for 
the modes in the previous section. However, now one 
has to consider hadronic matrix elements (J/ifrfo\Oi\B s ) 
with 0, four-quark operators, the calculation of which 
is a nontrivial task. In order to estimate the size of the 
B s — > /0-//1A decay rate, we use the generalized factor- 
ization approach, in which such quantities are replaced 
by products of matrix elements that are expressed in 
terms of meson decay constants and hadronic form fac- 
tors. The Wilson coefficients (or appropriate combina- 
tions of them) are regarded as effective parameters to 
be fixed from experiment. Using this ansatz, the decay 
amplitude of B s — > J/if/fo reads 

2G F l 

V2 

where e is the J /iff polarization vector, ps, the B s mo- 
mentum and fji$ = (416.3 ± 5.3) MeV the J/tf/ decay 
constant, ax is a combination of Wilson coefficients that 
can be extracted from <BR(B -> J/ifrK) H, assuming 
that ci2 is the same in the two processes. This requires 
the form factor F^ K . We use two different parame- 
terizations, obtained by short-distance (CDSS) 112011 and 
light-cone QCD sum rules (BZ) I2H1 . The result for 
the two sets of form factors is: \af DSS) \ = 0.394+g;^, 
| fl (#z)| _ q 25 + 0.03. We use the average value ax - 
0.32 ±0.11 and our result for the B s — > /o form factors 
to compute SKiBs — > J/ftfo), obtaining 



A(B S -> J/ifffo) = z -^-VcbV* s a 2 m^fj^Fi(mtJ(e*p Bs ) 



BK(B S -> J/tfff ) = (3.1 + 2.4) x 10" 



(12) 



large enough to be measured; notice that &R{B S — > 
J/ip(f>) = (1.3±0.4)xl0~ 3 [19]. Comparing these results 
to $¥R(B S — > JI^liPl) (L denotes a longitudinally polar- 
ized meson) computed using factorization, we find: 

R/u = ^ r^r = 0.13 ±0.06. (13) 

m !BK(B S ^ J/ifr L< f>L) 

A compatible result: — 0.2 - 0.3 was found in JfJ, 
using the ratio of D s decay widths to /o7r + and <p7i + . 

These considerations show that B s — > J/if/fo can be 
used to measure B s , since a large number of events is 
expected and it does not require an angular analysis to 
separate different CP components of the final state. This 
is also the case of B s — > Xa)4>, modulo the difficulty of 
the XcO reconstruction. Although suppressed in naive 
factorization, its branching fraction may be enhanced 
by non factorizable mechanisms i22ll as for B — > XcoK. 
On the basis of S U(3)f symmetry, we expect < B'R(B S — > 
XcQ<P) — O{\0 A ) as in the case of B — > XcoK* 112311 . 



5. Conclusions 

Exploiting the LCSR calculation of B s — > /o form fac- 
tors we find that the branching ratios of B s — > fo£ + f~ 
and B s — > fyw will be accessible at future machines, 
like a Super B factory, and at the LHCb experiment. 
We also predict SH(B S -> J/i//f )/S'R(B s J/ifi(f>) = 
0.13 + 0.06, thus B s — > J/if/fo is promising to access B s . 
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